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ABSTRACT 

Aims. We perform a quantitative morphological comparison between the hosts of active galactic nuclei (AGN) and quiescent galaxies 
at intermediate redshifts (z ss 0.7). The imaging data are taken from the large HST/ACS mosaics of the GEMS and STAGES surveys. 
Our main aim is to test whether nuclear activity at this cosmic epoch is triggered by major mergers. 

Methods. Using images of quiescent galaxies and stars, we created synthetic AGN images to investigate the impact of an optical 
nucleus on the morphological analysis of AGN hosts. Galaxy morphologies are parameterized using the asymmetry index A, the 
concentration index C, the Gini coefficient G and the M 10 index. A sample of ~ 200 synthetic AGN was matched to 21 real AGN in 
terms of redshift, host brightness, and host-to-nucleus ratio to ensure a reliable comparison between active and quiescent galaxies. 
Results. The optical nuclei strongly affect the morphological parameters of the underlying host galaxy. Taking these effects into 
account, we find that the morphologies of the AGN hosts are clearly distinct from galaxies undergoing violent gravitational interac- 
tions. Indeed, the host galaxy distributions in morphological descriptor space are more similar to undisturbed galaxies than to major 
mergers. 

Conclusions. Intermediate-luminosity (L x S 10 44 erg/s) AGN hosts at z ~ 0.7 show morphologies similar to the general population 
of massive galaxies with significant bulges at the same redshifts. If major mergers are the driver of nuclear activity at this epoch, the 
signatures of gravitational interactions fade rapidly before the optical AGN phase starts, making them undetectable on single-orbit 
HST images, at least with usual morphological descriptors. This could be investigated in future synthetic observations created from 
numerical simulations of galaxy-galaxy interactions. 

Key words, galaxies: active - galaxies: interactions - galaxies: evolution 



1. Introduction 

For a long time, active galactic nuclei (AGN) were considered a 
very special class of objects, characterized by an accreting super- 
massive black hole (SMBH) in their center. Only during the past 
decade, an SMBH has been found in practically all galaxies that 
are sufficiently close to determine the stellar kinematics in the 
center at high spatial resolution and that have a significant bulge 
component. It also became evident that the properties of SMBHs 
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and their host galaxies are coupled. The mass of the SMBH 
for instance is observed to be correlated with the luminosity 
(e.g. Marconi & Hunt 2003), velocity dispersion (e.g. Gebhardt 
et a l.l2000l Ferrarese & Merritt 120001 . and mass (e.g. Hiiring & 
Rix|2004) of the spheroidal component of its host galaxy. These 
correlations probably point towards a fundamental link between 
the formation and evolution of galaxies and their central black 
holes (e.g., Hopkins et al. [2008 ). However, some authors have 
demonstrated that the correlations might also arise from galaxy 
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mergers without any physically coupled growth of black holes 
and galaxies (e.g. Peng |2007l Jahnke & Maccio EoTTT) . 

The main mechanism that turns a quiescent SMBH into an 
AGN is still under debate. Candidates for this are secular pro- 
cesses such as bar-driven gas inflows, gravitational interactions, 
or galaxy mergers (for reviews, see Martini 2004] Jogee[2006). 
It is likely that the main triggering process will depend on the 
luminosity regime and/or the cosmic epoch under investigation. 

One possible approach for identifying this process is to com- 
pare the morphological properties of AGN host galaxies and 
quiescent (non-AGN) galaxies. To this end, some previous stud- 
ies have used descriptors such as CAS (Conselice et al. [2000). 
Grogin et al. (2005 ) for example found no difference between ac- 
tive and inactive galaxies at 0.4 < z < 1.3 in the Chandra Deep 
Fields North and South. Similarly, Pierce et al. (2007) used the 
Gini coefficient G and the M20 index (Lotz et al. 2004) to quan- 
tify the morphologies of X-ray selected AGN at 0.2 < z < 1 .2 in 
the DEEP2 survey and concluded that the majority of the AGN 
reside in nonpeculiar early-type hosts. These studies could be 
taken as evidence that signs of gravitational interactions in AGN 
at intermediate redshifts are scarce. However, it should be noted 
that these analyses did not include a decomposition of the AGN 
images into the nucleus and host components. Although an X- 
ray selection yields a high fraction of absorbed or so-called op- 
tically dull AGN, it is not clear whether optical nuclei had an 
impact on these results. 

In a different approach, some AGN studies focused on the 
environment. Serber et al. d2006) have found for instance that 
the galaxy density around quasars at z < 0.4 is enhanced with 
respect to inactive L* galaxies. The reported density excess is 
highest for the brightest quasars and separations R < 100 kpc. 
Hennawi et al. (2006) derived the autocorrelation function of 
quasars at 0.5 < z < 3 and showed an order-of-magnitude excess 
at R < 40 h~ l kpc in comparison to the expectation from large 
scales R > 3 h~ l Mpc. Myers et al. (120071 ) found a lower, albeit 
significant enhancement in the range 0.4 < z < 2.3. Silverman et 
al. (1201 11 1 compared galaxies at 0.3 < z < 1.1 with and without 
a companion galaxy at less than 75 kpc separation and showed 
an enhanced AGN fraction (by a factor of 2.6) in close pairs. 
According to a study by Ellison et al. (1201 11 1 based on > 10 s 
galaxies from the Sloan survey, this result also holds in the local 
universe. Although there is some disagreement on the magnitude 
of the density enhancement around AGN, these findings favor 
gravitational interactions or mergers as a triggering mechanism. 

An active nucleus with significant contributions to the to- 
tal flux will bias any morphological analysis of the host galaxy 
in the sense that it appears too bulge-dominated (e.g., Pierce 
et al. 1201 Oi l. Several studies have investigated ways to correct 
for the effect of a central point source (Sanchez et al. 120041 
Simmons & Urry [2008] Kim et al. 120081 Gabor et al. 120091 
Cisternas et al. 1201 11 1. Performing a decomposition into the host 
and nucleus components, Sanchez et al. (2004) found that 12 
out of 15 AGN at 0.5 < z < 1.1 reside in early-type galax- 
ies, with roughly one third of the hosts showing signs of grav- 
itational interactions. Moreover, the rest-frame U — V colors 
of the majority of the early-type hosts were significantly bluer 
than those of non-AGN early-types, indicating recent or ongo- 
ing star formation. The same has been found at similar redshifts 
for X-ray selected AGN in the Extended Chandra Deep Field 
South (Bohm et al. 120071 Silverman et al. |2008T >. Rafferty et 
al. d201 11 1 covered the wide redshift range 0.3 < z < 3.0 and 
showed evidence for a correlation between star formation rate 
and AGN fraction, which reaches 30 % for galaxies forming stars 
at if/ » 1000M o /yr. These findings could be understood in a sce- 



nario where some mechanism triggers star formation and, po- 
tentially with some delay, AGN activity, which eventually sup- 
presses star formation via feedback processes (e.g. Schawinski 
et al. |2007l Hopkins 12013 • 

Recent studies presented evidence against major mergers be- 
ing the main triggering process: Cisternas et al. d201 11 1 for in- 
stance found no strong distortions in > 85 % of their sample of 
AGN hosts at 0.3 < z < 1.0. Schawinski et al. (1201 It showed 
that, even at much earlier cosmic epochs up to z ~ 3, some 80 % 
of moderate-luminosity AGN (with Lx < 10 44 erg/s) reside in 
disk-like hosts, making a recent similar-mass merger very un- 
likely. 

In our analysis, we aim to further investigate the question 
whether AGN activity is connected to strong tidal interactions. 
Our sample covers intermediate redshifts (0.5 < z < 1.1) and X- 
ray luminosities (10 42 erg/s < Lx < 10 44 erg/s). We will for the 
first time use real galaxy images for the simulation of AGN im- 
ages and use these for a quantitative morphological comparison 
between quiescent and active galaxies at intermediate redshifts. 
At variance with previous studies, our sample is not X-ray but 
optically selected for AGN with broad emission lines (— » type- 
1). 

The paper is organized as follows: In Sect. [2] we introduce 
the surveys we used and Sect. [3] describes the selection of the 
AGN sample. The construction of the non-AGN comparison 
sample and the simulated AGN images is outlined in Sect. [4] 
We perform the morphological analysis of both samples and dis- 
cuss our results in Sect. [5j Sect. [6] will give a short summary. 
Throughout the paper, we assume a flat concordance cosmology 
with £2 A = 0.7, Cl m = 0.3 and H = 70kms _1 Mpc" 1 . Unless 
stated otherwise, all magnitudes are given in the AB system. 

2. Data sets 

Our morphological analysis relies on HST/ACS imaging data 
taken from the GE MS su rvey (Rix et al. EOOll l and the STAGES 
survey (Gray et al.|3p09 ). 

GEMS is a survey in the F606W and F850LP bands (similar 
to V and z), consisting of 78 ACS pointings. It covers a con- 
tinuous field of 28' x 28' in the extended Chandra Deep Field 
South. Details on the data reduction are given in Caldwell et 
al. ((2008 ). The STAGES survey extends over a slightly larger 
area (30' x 30') centered on the A901/2 multiple cluster system 
at z = 0.165. With 80 tiles of imaging through the F606W fil- 
ter, STAGES forms one of the largest mosaics ever taken with 
the HST Details on the data reduction can be found in Gray et 
al. ( 120091 . 

In brief, the data reduction comprised the standard steps 
(bias subtraction, flat-fielding, flux calibration), drizzling of 
the data from the original 0.05 arcsec/pixel to a final scale of 
0.03 arcsec/pixel, background estimation and variance determi- 
nation for each pixel. For both data sets, the individual exposures 
were multi-drizzled with a Gaussian kernelQ. For the purpose of 
this analysis, we used only the F606W images of GEMS, to al- 
low a combination with STAGES in a homogeneous data set. 

Both the STAGES and G EMS areas are covered by 
COMBO- 17 (Wolf et al. 120031 , a photometric redshift survey 
based on imaging in 12 medium- and 5 broad-band filters with 
the Wide Field Imager at the MPG/ESO 2.2m telescope on La 



1 The initial drizzling of the GEMS data was based on a square ker- 
nel. We here used data drizzled with a Gaussian kernel since it offers a 
significantly higher sensitivity in detecting faint host galaxies of AGN 
with brighter nuclei (Jahnke 2008, priv. com.). 
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Silla, Chile. COMBO- 17 provides the estimated redshifts and 
spectral energy distribution (SED) classifications of several 10 4 
objects. For galaxies brighter than R — 23 (Vega), the error on 
the photometric redshifts is typically cr z ~ 0.02(1 + z), (Wolf et 
al. I20081 I. For non-AGN galaxies, stellar mass estimates M„ are 
also available (Wolf et al. 120051 Borch et al . 120061 . 

3. AGN sample 

3. 1 . Local point spread function 

The analysis of the underlying host galaxy in AGN images is a 
challenging task, particularly in cases of low flux contrasts be- 
tween the host and nucleus. The key issue for accurately remov- 
ing the central point source and determining the host properties 
is a good knowledge of the point spread function (PSF). Its width 
and shape depend on the filter, the CCD position, the observation 
date, and the spectral shape of the source. Tests on GEMS data 
showed that the temporal variations are negligible compared to 
the spatial variations (see Sanchez et al. l2004t . 

Owing to the variations of the PSF shape across the CCD 
(with pixel-to-pixel flux variations of up to 20%), it is manda- 
tory to construct individual, local PSFs for each AGN. These 
local PSFs were determined by normalizing and averaging the 
nearest 35 stars from all tiles (only using stars from GEMS for 
GEMS AGN and from STAGES for STAGES AGN) around a 
given AGN position on the CCD. Very faint (V > 24) or sat- 
urated stars were rejected. The typical maximum distance be- 
tween the AGN and the stars that were included was ~ 40 arcsec 
for GEMS and ~ 25 arcsec for the STAGES data. 

A small fraction of the STAGES images (8 tiles out of 80) 
were obtained ~6 months later than the main block, which 
spreads over only three weeks. The tiles with delayed imaging 
(nos. 29, 46, 75 . . . 80; see Heymans et al. 2008 for details) show 
significant deviations of their PSFs from those of the data taken 
earlier and were therefore neglected in our analysis. To construct 
local PSFs just from these tiles was no option, it would have pro- 
duced PSFs with a much lower S /N. 

3.2. Definition of AGN sample 

Within the COMBO-17 catalog, there are 278 objects classified 
as AGN for which ACS imaging from STAGES or GEMS is 
available. We restricted our analysis to AGN in the redshift range 
0.5 < z < 1.1. The lower limit was chosen simply due to the lack 
of AGN identified with COMBO-17 at lower redshifts. The up- 
per limit is mainly motivated by the shift of the F606W filter's 
wavelength range into the rest-frame UV, which would strongly 
affect the morphological appearance and limit the comparability 
to objects at lower redshifts. Moreover, there is a dearth of quies- 
cent galaxies in COMBO- 1 7 at redshifts z > 1 . 1 usable as a com- 
parison sample. Because of the SED-based selection technique, 
the vast majority of the COMBO-17 AGN are type-1, i.e. they 
show broad emission lines and a prominent nucleus in optical 
images. We introduced a cut at R < 24 Vega mag for SED and 
redshift reliability. Finally, we rejected objects that were either 
located close to the edges of the tiles or observed during the de- 
layed STAGES observations (see previous section). This sample 
comprised 28 AGN with robust photometric redshifts and good 
ACS imaging data at 0.5 < z < 1.1. 

As a final constraint, we checked which of the AGN have a 
detectable host galaxy. To this end, the local PSF of each AGN 
was normalized to the AGN's flux within an annulus of a two- 
pixel radius and then subtracted from the AGN image. By defini- 



tion, this approach is an oversubtraction of the nucleus because 
the normalization to the central flux includes the contribution of 
the host galaxy. In other words, a scaled PSF subtraction pro- 
vides a conservative lower limit on the host-to-nucleus flux ratio 
H/N. Only objects showing a flux residual of more than 10% of 
the total flux were included in our final sample. This resulted in 
21 AGN with a median redshift of <z> = 0.71. 

3.3. Decomposition of AGN images 

We determined the H/N and the structural properties of the AGN 
host galaxies using the GALFIT package (Peng et al. [2002). 
GALFIT allows multi-component fitting of galaxy images in 
a x 1 minimization approach. To decompose the AGN images, 
the individually constructed local PSF was used as a model for 
the unresolved nucleus of the AGN, and a Sersic profile (Seisic 
1968) as a model for the host galaxy. This model can be used to 
approximate various types of surface brightness profiles with a 
Sersic index n — 1 being equivalent to an exponential disk and 
n — 4 to a de Vaucouleurs profile. 

The two-component GALFIT modeling introduces a total 
of ten free parameters. For the Sersic component, these were 
the position x s , y s , total brightness Vh os t, effective radius r e g, 
Sersic index n, position angle 6, and axis ratio q. For the nu- 
cleus model, only position x n , y n and brightness v nuc were free 
parameters. The user has to define initial guesses for all parame- 
ters, which GALFIT uses to start the fitting process. These initial 
guess values of m s , x s , y s , r e g, 0, m n , x n , y n were determined with 
Source Extractor (Berlin & Arnouts 1996). The flux ratio be- 
tween host and nucleus was set according to the estimate from 
scaled PSF subtraction. The initial guess value of the Sersic in- 
dex was n = 2.5, representing an intermediate type between an 
n — 1 exponential disk and a de Vaucouleurs profile. Fits were 
performed on thumbnails of size 256 x 256 pixel 2 , correspond- 
ing to 7.68 x 7.68 arcsec 2 . Neighboring objects and bad pixels 
were masked. 

GALFIT allows one to set constraint intervals on the free 
fitting parameters. We required the relative positions Ax = x s -x n 
and Ay = y s - y n to be within five pixels, constrained the Sersic 
index to 0.25 < n < 8 and the effective radius to 0.1 arcsec 
< r e ff < 3 arcsec. All fits were visually checked, in particular 
for the fit residuals. For a few objects, it was necessary to set 
additional constraints, e.g. on the position of the components, to 
acquire the best fitting results. 

We find median values of (v host ) = 23.05, (H/N) = 0.63, 
(n) = 1.51 and (r e s) - 0.22 arcsec corresponding to 1.6kpc at 
the average redshift (z) = 0.71 of the sample. Host parameter 
distributions and a test of their reliability are discussed in the 
following section. 

3.4. X-ray luminosities 

Our spectrophotometric classification rests on the COMBO-17 
data. By definition, all objects identified as AGN in this sur- 
vey (in the rest-frame UV/optical) are type-1, i.e. unabsorbed 
AGN. It is not clear whether our sample covers the same range 
in X-ray luminosities as X-ray-selected samples at similar red- 
shifts. Nuclear activity of different luminosity regimes might be 
triggered by different processes (e.g. secular evolution vs. grav- 
itational interactions), and this could complicate a comparison 
between our study and others. This section presents an estimate 
of the X-ray luminosities for our sample. 
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Ten objects in our sample are covered by the Chandra Deep 
Field South (CDFS) point source catalog (Lehmer et al. 2005 ). 
Their hardness ratios between the 0.5-2 keV and the 2-8 keV 
band (HR = (H - S)/(H + S), where H and S represent the hard 
and soft bands, respectively) cover the range -0.64 < HR < 
-0.40, confirming that they are unabsorbed. Adopting an X-ray 
power law slope of F = 2, we derived the X-ray luminosity at 
2-8 keV. They span values between Lx = 1.7 x 10 42 erg/s and 
L x = 3.2 x 10 44 erg/s with a median of (L x ) = 2.1 X 10 43 erg/s. 
For the other 1 1 objects, we had to apply an indirect estimate. 
Using the AGN template by Berk et al. (2001), we transformed 
the observed v-band nuclear brightnesses into the rest-frame B- 
band luminosities. These, in turn, were computed into X-ray lu- 
minosities following Hopkins et al. (2007 ). They cover the range 
5.9 x 10 41 erg/s < L x < 6.3 x 10 43 erg/s with a median of 
{Lx) — 7.2 x 10 42 erg/s. We hence find that, even though the 
sample consists of optically selected type-1 AGN with promi- 
nent nuclear sources, most of their X-ray luminosities are in the 
intermediate range. Only 3 out of 21 are powerful AGN with 
Lx > 10 44 erg/s. As a test, we used the approach of Hopkins 
et al. also to estimate the luminosities of the ten AGN covered 
by the CDFS. The indirectly computed X-ray luminosities differ 
from the CDFS-based ones by only ~ 0.22 dex on average. 

We now compare the luminosity range of our sample to 
X-ray-selected studies at similar cosmic epochs. Grogin et 
al. ( 120051 ) only gave a lower limit of Lx > 10 42 erg/s at the high- 
redshift end of their 0.4 < z < 1 .3 sample. The data sets of Pierce 
et al. (120071 0.2 < z < 1.2) Gabor et al. (120091 0.3 < z < 1.0), 
and Cisternas et al. (1201 II 0.3 < z < 1.0) span the ranges 
10 41 erg/s £ L x ^ 10 44 erg/s, 10 42 erg/s SL x S 4 x 10 44 erg/s 
and 2 x 10 42 erg/s ^ Lx ^ 3 x 10 44 erg/s, respectively. Despite 
of the differences in selection, our sample hence covers similar 
luminosities as these previous studies. 

4. Comparison sample 

In contrast to many previous studies we aim to take into account 
the impact of the optical nucleus on the quantitative AGN host 
galaxy properties. Thus we have constructed a special kind of 
comparison sample in three main steps. 

First, we defined a set of non-AGN galaxies that is matched 
to the redshifts and magnitudes of our sample of 21 AGN hosts. 
Secondly, we created synthetic AGN images from these quies- 
cent galaxies that cover the same H/N as the real AGN (sec- 
tion |43]). Finally, the synthetic AGN images were decomposed 
into the host and nucleus components using GALFIT in the same 
way as the real AGN. After subtracting the model nucleus from 
the AGN image, the recovered host galaxies of the synthetic 
AGN were combined with the host galaxies of the real AGN 
and analyzed in concentration, asymmetry, Gini and M20 index 
(see Section|5]l. 

4.1. Selection 

The basis for the selection of the comparison sample was the 
COMBO- 17 catalog of the GEMS field. For the pre-selection, 
the following constraints were used: 

1. Apparent Vega magnitude R < 24. In this magnitude regime, 
we can rely on photometric redshifts accurate to within 3- 
4% a nd rob ust SED classifications from COMBO- 17 (Wolf 
et al. 120031 . 

2. Source classified as a "galaxy" on the basis of its SED (but 
not as an AGN). 
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Fig. 1. Comparison between the redshifts of the AGN (dashed) 
and the sample of quiescent galaxies (solid). The histograms are 
normalized to the same total number of objects. 

3. No X-ray counterpart in the point source catalog of the 
Extended Chandra Deep Field South (Lehmer et al. 2005 ). 
This criterion was applied to exclude absorbed or "optically 
dull" AGN. 

4. Sersic index n > 1.5 (derived on the F850LP images, i.e. z- 
band, by Haussler et al. 2007). This constraint was set to 
restrict the sample to galaxies with significant bulge com- 
ponents. According to the relation between bulge mass and 
SMBH mass, host galaxies with small bulges or even without 
a bulge are not expected to be present in the AGN sample. 

These criteria were fulfilled by 3536 galaxies with GEMS 
data. In the next step, ten quiescent counterparts were selected 
for each of the 21 AGN hosts within a bin Az/z < 0.025 of the 
AGN redshift and within Av host < 0.05 mag of the total bright- 
ness of the host as determined with the GALFIT decomposition. 
When less than 10 non-AGN were found in the catalog within 
these limits, the limits were slightly relaxed and the search re- 
peated iteratively until ten counterparts could be selected. After 
visual inspection of all 210 objects in the resulting sample, seven 
galaxies had to be rejected because they were unresolved. The fi- 
nal non-AGN sample hence contained 203 galaxies. 

The COMBO photometric redshift distributions of the 21 
AGN hosts and the 203 quiescent galaxies are shown in Fig.Q] 
The z-distribution of the AGN was normalized to a total of 
203 objects to allow a direct comparison between the two sam- 
ples. The median photometric redshift of the AGN sample is 
(z) = 0.71, virtually identical to the (z) = 0.70 of the quiescent 
sample. Fig.|2]displays the distributions in total apparent V mag- 
nitude of the AGN hosts (from GALFIT) and the quiescent com- 
parison sample (Source Extractor MAG_AUTO values). Again, 
both samples' distributions are normalized to the same number 
of objects. The median total magnitudes of the two samples are 
(Vhost) = 23.05 (AGN) and (vh os t) = 23.33 (comparison sample). 
The slight difference can be explained with the relaxation of the 
constraints partly necessary during the selection of the compari- 
son sample. 

4.2. Visual classification 

To calibrate the following quantitative morphological analysis, 
we performed a visual classification of all galaxies in the qui- 
escent sample. We decided to use a very plain scheme. Its 
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z=051 undisturbed z=0.52 undisturbed z=0.G7 undisturbed z=0.70 undisturbed z=0.86 undisturbed z=0.87 undisturbed 





z=0.57 unclear z=0.S8 unclear z=0.34 unclear z=0.95 unclear z=1.03 unclear z=1.03 unclear 



Fig. 3. HST/ACS v-band images of galaxies from the non-AGN comparison sample. Each row displays six examples of the three 
different morphological types that were denoted via visual classification (with at least four out of eight classifiers agreeing). The 
bar in each image has a length of one arcsec, corresponding to ~ 6 kpc at z — 0.5 or ~ 8 kpc at z = 1 . 
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Fig. 2. Distribution of the AGN host galaxy magnitudes de- 
rived with GALFIT (dashed) compared to the quiescent sample 
(solid). The histograms are normalized to the same total number 
of objects. 



main purpose is to test whether major mergers are the triggering 
mechanism of AGN at redshifts 0.5 < z < 1 . 1 . To this end, we 
distinguished between three morphological classes: 

1. Major mergers: galaxies that are heavily distorted ("train- 
wrecks"), show two nuclei of similar luminosity, or two pro- 
nounced tidal tails of equal length. 

2. Undisturbed galaxies: no signs of recent or ongoing grav- 
itational interactions. 

3. Unclear cases: all other objects, including minor merg- 
ers, ambiguous major merger candidates, or morphologies on the 
edge between undisturbed and mildly disturbed. 



Using this simple scheme, all 203 galaxies in the compar- 
ison sample were visually inspected and classified by eight of 
the authors (AB, DB, MaB, MG, GH, SJ, DM, AR). Comparing 
the individual classifications, agreement was found between four 
classifiers or more in 88% of the cases; between five or more in 
84% of the cases; between six or more in 67% of the cases, and 
between seven or more in 48% of the cases. 

The classifications were merged with a constraint of at least 
four classifiers agreeing on a certain type. Cases with less than 
four identical classifications were considered unclear. If an ob- 
ject was classified as a major merger four times and as unclear 
four times, it was also defined as unclear. We refer to the crite- 
rion of at least four identical classifications per object as the 4/8 
sample. Examples are presented in Fig. [3] Similarly, a scheme 
5/8 requires that five or more classifiers agree (or else, the given 
object is unclear); and equivalently for the 6/8 and 7/8 samples. 
The respective number of major mergers, undisturbed, and un- 
clear objects are summarized in TableQ] 

We tested how well the various combined classifications dis- 
tinguish between undisturbed galaxies and major mergers on the 
basis of morphological descriptors such as concentration C and 
asymmetry A, to choose a scheme for our final analysis. This is 
described in Sect.P 



4.3. Synthetic AGN images 

The modifications of the quiescent galaxy ACS images de- 
scribed in the following constitute the backbone of our analy- 
sis. The same kind of residuals that occur after subtracting the 
nucleus in the central part of an AGN host were introduced in 
the quiescent galaxy images. This ensures a direct comparabil- 
ity of the morphologies of AGN hosts and non-AGN galaxies. 
We achieved this by adding a synthetic nucleus to a given quies- 
cent galaxy, constructing its local PSF and finally decomposing 
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Table 1. Visual classifications of the 203 quiescent comparison 
galaxies. 



Criterion 


Maj. Merg. 


Undist. 


Unclear 


(1) 


(2) 


(3) 


(4) 


4/8 


15 


95 


93 


5/8 


14 


95 


94 


6/8 


8 


86 


109 


7/8 


4 


68 


131 



Note: Col. (1) denotes the required minimum number of classifiers 
agreeing on the type of a given object. Cols. (2) to (4) compare the 
respective numbers of major mergers, undisturbed and unclear objects. 



it into the host and nucleus component with GALFIT. Thus, we 
analyzed synthetic and real AGN in exactly the same way. For 
the synthetic nuclei, we used real stellar images from GEMS. 
It is hence ensured that the deviations between the light profile 
shape of a nucleus and the corresponding local PSF will be sta- 
tistically the same as in the case of the real AGN, and the resid- 
uals of the nucleus in simulated AGN and real AGN images — 
after subtracting the best-fit model nucleus — will have the same 
noise properties. 

The stars we used as synthetic nuclei were selected from the 
GEMS catalog, constrained to the v-band magnitude range of 
the nuclei of the real AGN, 19.6 < v < 25.3. Since the PSF 
shape is sensitive to the SED of a source and the nuclei of AGN 
have very blue colors, we set color cuts of (B — V) < 0.6 and 
(V — R) < 0.6 (both Vega magnitudes). Five stars were rejected 
because they showed cosmics or bad pixels close to the central 
peak. The final list for creating synthetic nuclei consisted of 58 
stars with median colors of (B -V) ~ 0.44 and (V - R) ~ 0.43 
(Vega magnitudes). 

For each of the quiescent galaxies, a star was selected from 
this list such that the flux ratio between galaxy and star was the 
same as the H/N of the AGN to which the quiescent galaxy ini- 
tially was assigned in the construction of the comparison sam- 
ple. Small normalizations of about a few percent were applied 
to the stellar fluxes to achieve exactly the same flux ratios be- 
tween a galaxy and the star on the one hand and AGN host and 
nucleus on the other. After normalization, the quiescent galaxy 
(synthetic host) and stellar image (synthetic nucleus) were sim- 
ply co-added with the star positioned at the galaxy's central coor- 
dinates according to the GALFIT modeling on the initial image. 
The variance frames of the two components were also combined 
to use them in the subsequent GALFIT decomposition. 

The local PSF for each synthetic AGN was created based on 
the CCD position of the star that was used as a synthetic nu- 
cleus. The 35 nearest stars around the CCD position of a given 
synthetic nucleus were normalized and averaged. Local PSFs for 
the synthetic AGN hence have the same quality as those of the 
real AGN. Finally, using GALFIT, all synthetic AGN were de- 
composed into nucleus and Sersic host with the same configura- 
tion as described for the real AGN (section l3~3l l. The main steps 
of the creation and decomposition of the synthetic AGN are il- 
lustrated in Fig. [4] 

The H/N values of the real AGN and the synthetic AGN 
are compared to each other in Fig. [5] As in the case of the red- 
shifts and host magnitudes, the two samples differ only slightly 
in H/N, with median values of (H/N) = 0.57 and (H/N) = 0.63 



♦ • 




Fig. 4. Illustration of the AGN simulation procedure. Upper left. 
original quiescent galaxy at redshift z = 0.52, apparent bright- 
ness v = 23.27. Upper middle: star of brightness v = 21.46 used 
as a synthetic nucleus. Upper right: simulated AGN (co-added 
galaxy and star) with a host-to-nucleus ratio of H/N = 0.18. 
Lower left: GALFIT model of the host galaxy (from fitting a 
nucleus + host model to the simulated AGN). Lower middle: re- 
covered host galaxy after subtracting the GALFIT model of the 
nucleus from the simulated AGN image. Lower right: residual 
after subtracting the combined galaxy + nucleus model from the 
simulated AGN image. Each image is ~ 2 arcsec on the side. 



80 - 



40 



H/N 

Fig. 5. Distribution of the host-to-nucleus ratios H/N of the 
real AGN (dashed) compared to the simulated AGN (solid). The 
histograms are normalized to the same total number of objects. 



for the simulated AGN and real AGN, respectively. These devi- 
ations do not stem from the selection process itself but from the 
limited accuracy of the GALFIT decomposition; we discuss this 
in detail in the next section. 

The visual classifications described in the previous section 
were also applied to the real AGN sample after subtracting the 
best-fit models of the nuclear point sources from the images. 
These classifications were made as a blind test with all nucleus- 
subtracted real AGN and all nucleus-subtracted simulated AGN 
randomly mixed. We present an overview in Table [2] It is strik- 
ing that none of the real AGN is classified as a major merger. 
This is discussed in more detail in Sect. 15.41 in combination with 
the results from our quantitative morphological analysis. All real 
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Fig. 6. Visual classifications of the nucleus-subtracted real AGN sample (with subtracted model nucleus from two-component 
modeling as described in Sect l3.3b . The bar in each image has a length of one arcsec. 



Table 2. Visual classifications of the 21 real AGN hosts. 



Criterion 


Maj. Merg. 


Undist. 


Unclear 


(1) 


(2) 


(3) 


(4) 


4/8 





10 


11 


5/8 





9 


12 


6/8 





7 


14 


7/8 





3 


18 



Note: Columns are the same as in Table[T] 



AGN and their respective visual classes (based on the 4/8 crite- 
rion) are shown in Fig|6] 

4.4. Stability of the AGN decomposition 

Our synthetic AGN images can be used to test the GALFIT re- 
covery of the main host galaxy parameters such as total bright- 
ness, effective radius, and Sersic index for a broad range of H/N. 
Similar tests have been carried out by Sanchez et al. (2004), 
Kim et al. (120081 1. and Gabor et al. (120091 1: but in these studies, 
no real galaxy images but simulated, perfectly symmetric hosts 
have been used. Only Simmons & Urry (2008 ) used observed 
data for part of their tests. 



Even though studying the robustness of AGN host galaxy 
parameters derived with GALFIT is not the raison d'etre of our 
analysis, it is an interesting byproduct. The plots in Figs. Q and 
|8]compare the initial host parameters (with suffix ini), as derived 
on the images prior to adding a synthetic nucleus, and the re- 
covered host parameters (suffix sim) from simultaneous fitting 
of galaxy and synthetic nucleus. 

In Fig. [7^, we show the total brightness Vh os t,ini of the qui- 
escent galaxies (single Sersic modeling on the original images, 
before the addition of synthetic nuclei) plotted against the de- 
viations Avhost = Vhost.ini - Vhost.sim between the initial bright- 
ness Vhost,ini and the brightness Vhost.sim determined from two- 
component GALFIT modeling after adding the synthetic nu- 
cleus. The values of Avh os t are basically distributed around zero 
with a scatter of ~ 0.25 mag. The scatter is mostly indepen- 
dent of the galaxy magnitude except for objects brighter than 
Vhost.ini ~ 21.3. It is surprising that most of these galaxies have 
high positive values of Avh os t, meaning that the fluxes of these 
objects are overestimated when they are fitted simultaneously 
with a synthetic nucleus. It is quite obvious that this occurs be- 
cause the fitting algorithm attributes part of the nucleus flux to 
the host, because the Sersic indices of most of these galaxies 
are overestimated, too — this leads to low values of the frac- 
tion nj n i/n s i m between Sersic index in the original image and the 
simulated AGN, depicted by filled green circles; see also panel 
e). 

Fig. |7j3 shows the deviations between original host magni- 
tude and recovered host magnitude as a function of input H/N^. 
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Fig. 7. GALFIT recovery of host galaxy parameters for the simulated AGN sample, a) Deviations between input host magnitude 
Vhost,ini and recovered host magnitude as a function of input host magnitude, b) Deviations between original host magnitude and 
recovered host magnitude Vh os t,sim as a function of input H/N ia i. c) Deviations between input effective radius r e ff,i n i and recovered 
effective radius r e ff s i m as a function of input effective radius, d) Deviations between input effective radius and recovered effective 
radius as a function of input H/N. e) Deviations between input Sersic index n; n ; and recovered Sersic index n s j m as a function of 
input host magnitude./) Deviations between input Sersic index n; n j and recovered Sersic index n s ; m as a function of input H/N. 



The scatter in Avhost strongly increases toward low values of 
H/Ni n i, as expected. Moreover, in the regime H/N- m i < 1, i.e. in 
cases where the nucleus outshines the host, the recovered host 
magnitudes Vh 0s t,sim tend to be too bright. This trend indicates 
that part of the bright nucleus flux is attributed to the host com- 
ponent by the fitting algorithm. A similar results has been found 



by Kim et al. (2008). These authors argued that mismatches be- 
tween the nucleus profile and the input local PSF cause this sys- 
tematic effect. 

Fig. [7J: shows the input effective radii r e ff,mi determined on 
the original images plotted against the deviations between in- 
put effective radius r e ff,i n i and recovered effective radius r e ff jS im- 
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This parameter is accurate to within less than 30% for 3/4 of the 
objects. Fig. [7Jl shows the deviations between input and recov- 
ered effective radius as a function of input H/Nm- Like the host 
brightnesses, the effective radii have a larger scatter toward low 
H/N. 

In Fig. [7£, the ratio An = log(ni n ;/n s i m ) between original 
and recovered Sersic index is plotted as a function of input to- 
tal host magnitude Vh 0s t,ini- Evidently, n is a parameter that is 
much harder to determine than the brightness or effective ra- 
dius of the host galaxy. Toward fainter host magnitudes, the 
Sersic indices are systematically underestimated in a decom- 
position into nucleus and host component. For the distribution 
in H/N considered here (the median for our AGN sample is 
(H/N) = 0.63), this effect starts to kick in at host magnitudes 
fainter than Vh os t,ini ~ 23. Owing to this systematic effect and the 
very large scatter in An, it is impossible to use the best-fit Sersic 
index to determine the morphological type of an AGN host. A 
similar trend has been found e.g. by Sanchez et al. (f2p04) on the 
basis of simulated host galaxies and nuclei. Those simulations 
used symmetric single-component host galaxies (Sersic profiles 
with either n = 1 or n = 4). In contrast to this, many of the 
host galaxies in our sample consist of two or more components 
(bulge, disk, potentially spiral structure, etc.). This higher com- 
plexity of the host profiles leads to a much larger uncertainty 
in the recovered Sersic indices than found by Sanchez et al. A 
comparison to the results of Simmons & Urry (2008) is difficult 
because the authors did not use constraints on the Sersic index 
in their fits. 

Fig- Ef displays the ratio An = log(n; n j/n s j m ) between origi- 
nal and recovered Sersic index as a function of the input H/N. 
At low H/N, the scatter in An is approximately 0.4 dex. 

Fig. [8^ shows the host magnitude deviations Avhost plotted 
as a function of the Sersic index deviations An. Not surprisingly, 
there is an anti-correlation of these two parameters. A too faint 
model profile for the host galaxy (Avh ost < 0) is accompanied by 
an underestimated Sersic index, because a fraction of the host 
galaxies' flux is assigned to the model nucleus. Since this effect 
occurs at small galactocentric radii, the host model profile be- 
comes too shallow (vice versa for an overestimated host flux). A 
false host model flux results in a wrong value of the model H/N 
(see Fig. [SJ?, displaying AH/N as a function of the host magni- 
tude deviations Avhost)- Fig. [SJ: shows the ratio between the input 
H/N m i and the recovered value H/N s \ m , plotted as a function of 
H/N m i. The H/N is slightly underestimated for low input values 
and slightly overestimated for high input values, which partly 
is a consequence of the systematic error on Vh os t at low H/N. In 
Fig.[8jl, the ratio between the input and recovered H/N is plotted 
against the original total host magnitude Vh 0s t,ini- In general, H/N 
is overestimated for very bright hosts — which is induced by the 
frequently overestimated host fluxes in this regime (cf. Fig.|7h)- 

The comparisons between single-component light profile 
modeling of quiescent galaxies and two-component modeling of 
these galaxies with added synthetic nuclei provides a system- 
atic error in the decomposition of AGN images using GALFIT. 
Ignoring extreme outliers and taking only the best 95% of the 
fits into account, host fluxes, H/N, effective radii and Sersic in- 
dices are recovered to within fractional errors of cr v h Os t=0.23, 
c"h/n=0.34, and o" le ff=0.50, cr n =0.90, respectively. Since the 
scatter in host magnitude, effective radius and Sersic index is 
dominated by objects with low H/N (cf. Fig [T]}, d and f), we 
repeated this test by rejecting objects with H/N < 0.4. We then 
find respective fractional errors of cr vhost =0.18, ct H /n=0.46, and 
c re ff=0.22, cr n =0.82. In particular for the effective radius, these 
values are more representative for an AGN decomposition with 



not too faint hosts. The numbers also show that the Sersic indices 
of distant type-1 AGN hosts are very hard to recover. 

5. Morphological analysis and discussion 

Prior to the following analysis, the best-fit model nuclei from the 
GALFIT decomposition were subtracted from all synthetic and 
real AGN images. 

We define a galaxy's outer contours using the segmentation 
maps from Source Extractor because they determine which pix- 
els of an image are assigned to a given galaxy. All these maps 
were visually inspected. In a few cases, e.g. when the deblending 
of multiple galaxies was not satisfying, Source Extractor was re- 
run with different input parameters until the segmentation maps 
were accurate. In our approach, major merger candidates were 
not deblended into the individual sub-components but analyzed 
as one single object. In the SExtractor configuration used here, 
the flux at the outermost isophotes of a given galaxy corresponds 
to the 1 <t sky noise level. 

5. 1 . Definition of morphological descriptors 

We here rely on two parameter spaces that have been frequently 
used in the literature, in particular for distant galaxies. These are 
the concentration index C and asymmetry index A (Conselice 
et al. I20001 I. the Gini coefficient G and the M20 index (Lotz et 
al. l2004l l. 

The asymmetry parameter A is measured by subtracting a 
galaxy image rotated by 180° from the original image: 

, . /2 W |/-/i8ol\ . /Z W |B--BimI\ 

A = mm — -mm — ; . (1) 

Here, / is the original image, /iso is the image rotated by 180° 
about the galaxy's center, and B and B\^) are the background and 
rotated background. This latter contribution to A due to the back- 
ground noise is estimated from blank sky regions in the vicinity 
of the object under investigation. The sums are computed over all 
pixels within the lcr isophotes of the galaxy; B covers the same 
number of pixels as /. The minimization is computed by shifting 
the galaxy image on an integer pixel grid, with a maximum shift 
of eight pixels for our data sets. 

There are various definitions of the concentration index C in 
the literature that only have slight differences to each other. We 
adopted C as the ratio between two fluxes in elliptical areas G\ 
and G2, which have the same axis ratio and position angle but a 
different major axis with 02 = 0.3 a\ (Abraham et al. 1 1996k 

C=|^. (2) 

Here, /y is the flux at the pixel (i, j) within the inner isophotal 
region G2 and /# is the flux at the pixel (k, I) within the larger 
isophotal region G\ . 

The Gini coefficient has originally been developed for appli- 
cations in social sciences. When computed for galaxy images, 
it is a measure of the concentration of an object's light profile, 
similar to C, but insensitive to the definition of a central pixel. 
We adopted the Gini coefficient as given in Lotz et al. (2004): 

G= mr! n I>- w - Wil. (3) 
\ f\n(n -l)*f 

where \fi\ are the absolute flux values of the z'-th pixel (ranked in 
increasing order) and \ f\ is the mean of the absolute flux values 
of all n pixels assigned to an object. 
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Fig. 8. GALFIT recovery of host galaxy parameters for the simulated AGN sample, a) Deviations between input host magnitude 
Vhost,ini and recovered host magnitude Vh 0s t,sim as a function of the deviations between input Sersic index n; n j and recovered Seisic 
index n s j m . b) Deviations between input H/N ini and the recovered H/N sim as a function of the deviations between original host 
magnitude and recovered host magnitude, c) Deviations between input H/N ini and the recovered H/N sim as a function of input H/N. 
d) Deviations between input H/N ini and the recovered H/N sim as a function of input host magnitude. 



The M20 index is based on the total second-order moment 
M tot , which is the flux in each pixel /; multiplied by the squared 
distance to the center of the galaxy, summed over all pixels of a 
given object: 

n n 

M tot =Y i M i = Y i f i - {{ Xi - x c f + (y t - y c f), (4) 

i i 

where x c ,y c is the galaxy's center. The center is computed by 
finding x c ,y c such that M tot is minimized. M20 then is defined 
as the normalized second-order moment of the brightest 20% of 
the galaxy's flux. To compute M20, the pixels are rank-ordered 
by flux, and M, is summed over the brightest pixels until the sum 
of the brightest pixels amounts to 20% of the total galaxy flux, 
followed by a normalization to M tot (Lotz et al.|2p04): 

M 20 = log (j^j while J] fi < 0.2f tot , (5) 

Here, f tot is the total flux of a galaxy and /; are the fluxes for 
each pixel i, ordered according to their flux values such that fa 
is the brightest pixel, fa the second brightest, etc. Basically, M20 
traces the spatial distribution of bright features such as off-center 
star-forming regions. 



5.2. Quantitative morphological comparison of AGN hosts 
and quiescent galaxies 

Our visual classifications of the comparison sample will be 
used as a calibration in the C/A and G/M20 parameter spaces. 
Since these classifications were carried out prior to including a 
synthetic optical nucleus, they are much more robust than the 
nucleus-subtracted host galaxy images of the real AGN that are 
affected by the nucleus residuals in the center of the images (cf. 
Fig.0. 

First we defined when to call a visual classification reliable. 
Eight team members visually inspected the galaxies in the non- 
AGN comparison sample and classified them as a major merger, 
undisturbed, or unclear. Combining the various classifications 
for each object, we constructed various samples (as described in 
Sect. 14.2b . Only using cases with at least four identical classifi- 
cations, we obtained a sample 4/8, five or more identical clas- 
sifications define a sample 5/8, etc. We used a two-dimensional 
Kolmogorov-Smirnov test to quantify the distinguishability be- 
tween major merger candidates and undisturbed galaxies in the 
parameter spaces C/A and G/M20. Table [3] shows the respective 
probabilities P that both samples are drawn from the same par- 
ent distribution. The probability P increases with the minimum 
number of agreeing classifications. This mainly is an effect of the 
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Table 3. Two-dimensional KS test comparing major mergers and 
undisturbed galaxies in the quiescent comparison sample. 



Sample 


C/A 


G/M20 


(1) 


(2) 


(3) 


4/8 


1(T 6 


8 x i(r 6 


5/8 


2 x 1(T 6 


1.8 x lfr 5 


6/8 


4 x 1(T 5 


3.1 x 10 -5 


7/8 


i.i x itr 3 


9.1 x lO -4 



Note: Col. (1) denotes the required minimum number of classifiers 
agreeing on the type of a given object; Col. (2): probability P that both 
samples stem from the same parent distribution in C/A space; Col. (3): 
same for G/M20 space. 

Table 4. Two-dimensional KS test comparing major mergers 
and undisturbed galaxies in the simulated AGN sample (after 
removal of the central point sources). 



Sample 


C/A 


G/M 20 


(1) 


(2) 


(3) 


4/8 


8.6 x itr 4 


1.1 x lfr 3 


5/8 


5.7 x itr 4 


2.7 x icr 3 


6/8 


5.2 x itr 3 


3.9 x itr 2 


7/8 


1.1 x irr 2 


4.4 x icr 3 



Note: Columns are the same as in Table[3] 

sample sizes. As a test, we restricted all samples to the size of the 
7/8 scheme (by randomly selecting subsets from the other larger 
samples) and recomputed the KS test. The probabilities then are 
comparable between the four different schemes. However, since 
a larger sample offers more robust statistics, we decided to use 
scheme 4/8 in the following. 

Our main morphological analysis is based on galaxy im- 
ages with removed central point sources (cf. Fig. 0). For this 
reason we tested the distinguishability between major mergers 
and undisturbed galaxies also with the images of the nucleus- 
subtracted simulated AGN. As shown in Table |4] the contrast 
between the two morphologies is degraded with respect to the 
original images, because of the residuals in the galaxy centers 
after the point-source subtraction (just for completeness, we give 
the KS test results for all four classification schemes). 

In Fig. [9] (left-hand plot), we show the distribution of the 
quiescent galaxies in terms of C vs. A. Major merger candidates 
have systematically higher asymmetry indices A and lower con- 
centration C than undisturbed galaxies. This is well-known from 
observations as well as from synthetic images created from nu- 
merical simulations (e.g. Lotz et al. 2008). The galaxies with un- 
clear morphology are located between major mergers and undis- 
turbed galaxies, showing substantial overlap with the former. A 
number of galaxies that could not be clearly characterized visu- 
ally might thus be major mergers. Still, the KS test probability 
that major mergers and unclear objects stem from the same par- 
ent distribution is low: P = 0.004. 

In G/M20 space (right-hand plot of Fig. [9]), the situation is 
similar: major mergers populate a region of high M20 values (i.e. 



bright off-center features) and lower G indices compared to the 
undisturbed galaxies. The distinction between the two classes 
seems slightly worse than in C/A space; this indeed is reflected 
in the higher KS test P values for G/M^o with respect to C/A 
(Table©. 

The distributions in C/A and G/M20 are completely differ- 
ent in the presence of an optical nucleus. In Fig. [10] we show 
simulated AGN (quiescent galaxies plus synthetic nucleus) and 
real AGN (including their optical nuclei). Note that the displayed 
range in M20 is slightly different from Fig. [9] Compared to the 
original quiescent galaxies, the synthetic AGN have much higher 
concentration indices and Gini coefficients, which is an obvi- 
ous effect of the nuclear point source. The M20 indices become 
lower, indicating the loss of structure because the combined light 
profile of host+nucleus is smoother than that of the host alone. 
The mean asymmetry of all types (median value (A) s ; m = 0.3) 
is almost unchanged with respect to the original galaxy images 
( Worig = 0.27). This is somewhat surprising because the nuclear 
point sources should result in a less elliptical combined light 
profile. However, the synthetic nucleus is usually not perfectly 
centered on the pixel grid; just like the nucleus in an observed 
real AGN. Since the asymmetry index is minimized on an inte- 
ger pixel shift basis, even an almost symmetric light profile can 
therefore yield an untypically high value for A. One other net ef- 
fect of the nuclei was that the distribution in A became narrower. 
The distinguishability between the morphological classes of the 
simulated AGN hosts is, when compared to the initial distribu- 
tions in Fig. [9] clearly degraded. Repeating the KS test to quan- 
tify how well undisturbed and major mergers are separated, we 
find P = 0.001 in C/A and P = 0.004 in G/M 20 . Simulated and 
real AGN cover similar regions in the parameter spaces. Before 
we perform a quantitative comparison, we go one step further. 

For the next computation of C, A, G, and M20, the best-fit 
nuclei from two-component modeling with GALFIT were sub- 
tracted from the simulated and real AGN. The resulting distri- 
butions are shown in Fig.Qj] In C/A space, mergers and undis- 
turbed objects are better separated than before subtracting the 
nucleus (cf. Fig. [TOt . However, there clearly is a stronger over- 
lap between these two types than there was based on the original 
images (cf. Fig. [9j>. This is reflected in an increase of the KS test 
probability when comparing merging and undisturbed galaxies: 
values increase from P = 10~ 6 to P « 10~ 3 in C/A and from 
P * 10~ 5 to P * 10 3 in G/M 20 , as given in Tables and g] 
Some objects feature asymmetry indices above unity, which in 
the case of normal galaxy images would be impossible. This ef- 
fect indicates additional noise in the center of the images that is 
not related to the sky background but residuals from the nucleus 
subtraction. As we discuss below, this particularly affects objects 
with low H/N. 

We will now use the KS test to determine which subsam- 
ple of the quiescent galaxies is most similar to the AGN hosts. 
Or, more precisely, which subsample of the nucleus-subtracted 
synthetic AGN hosts is most similar to the real AGN hosts. 
Comparing the AGN to major mergers, undisturbed, and un- 
clear objects in C/A space, we find P mm = 0.015, P U ndis = 0.32 
and P unc i = 0.034, respectively. The distribution of the AGN 
hosts in asymmetry and concentration is by far best matched by 
the undisturbed galaxies, while unclear cases show the second 
highest — albeit much lower — probability. The same compar- 
ison in G/M20 space yields P mm = 0.040, P U ndis = 0.22 and 
''unci = 0.033, respectively. Again, the distribution of undis- 
turbed galaxies is most similar to the AGN hosts, while major 
mergers are much less likely to stem from the same parent dis- 
tribution. 
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Fig. 9. Distribution of the quiescent galaxies in concentration index C versus asymmetry A (left) and in Gini coefficient G versus 
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Fig. 10. Distribution of the simulated AGN (i.e. a quiescent galaxy with an additional synthetic nucleus (before subtraction), 
depicted here as circles) and real AGN host galaxies (open diamonds), in concentration index C versus asymmetry A (left) and in 
Gini coefficient G versus M20 (right). Tiny open symbols depict ambiguous cases, large open symbols show undisturbed galaxies, 
large filled symbols are classified as major mergers. With respect to the original galaxy images (the distributions of which are shown 
in Fig. [9j>, the distinguishability between different morphological types is degraded. 



We also compared the AGN hosts to the whole quiescent 
sample. The KS test then yields P. d \\ = 0.20 in C/A or P^i = 
0.12 in G/M20. In both parameter spaces, these probabilities are 
slightly lower than those comparing between AGN hosts and 
undisturbed galaxies. This indicates that the AGN hosts com- 
prise a mix of morphological classes that is slightly different 
from the comparison sample, in the sense that the former con- 
tain a higher fraction of undisturbed galaxies than the latter. 

These numbers argue against a significant fraction of ma- 
jor mergers among the AGN hosts. Frequent minor mergers, on 
the other hand, cannot be ruled out. The unclear category is 
a mixed bag, containing minor mergers as well as ambiguous 
cases that might partly be major mergers. If major mergers are 
not the main triggering mechanism, these ambiguous cases will 
result in a relatively low KS test probability, just like we have 
found. We tested an alternative visual classification scheme that 
distinguishes between minor and major mergers. To do this on 
an object-by-object basis is of course a challenging task. The us- 



age of this alternate scheme confirmed our results, albeit at lower 
statistical significance, on the basis of KS statistics. This is why 
we here adopted a simple scheme with only three different mor- 
phological types, as a straightforward test of the major merger 
scenario. 

The presented distributions in the parameter spaces spanned 
by C, A, G, and M20 already gave some insight into the effects 
of an optical nucleus in the morphological analysis. We now in- 
vestigate this in more detail for each of the four descriptors in 
Fig. [12] In the upper left plot, C m \ is the initial concentration 
index of a non-AGN galaxy before adding a synthetic nucleus, 
C s im gives the index after addition and subtraction of a synthetic 
nucleus. The difference between the two is normalized to the 
absolute initial value via AC = (C; n j - C s i m /|C m i|). This param- 
eterization allows us to distinguish between systematic trends, 
showing up as changes in the average value of AC, and a loss of 
information leading to an increased scatter in AC. We show AC 
as a function of H/N. 
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Fig. 11. Distribution of the simulated and real AGN host galaxies after subtracting the best-fit model nucleus in concentration index 
C versus asymmetry A (left) and in Gini coefficient G versus M20 (right). Tiny open symbols depict ambiguous cases, large open 
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are better separated than before subtracting the nuclei (Fig. ITOt. 



After subtracting the synthetic nucleus, the concentration in- 
dices C s i m are on average slightly underestimated for high H/N 
values, corresponding to AC > 0. This trend is the equivalent 
of the underestimated Sersic indices « s j m in Fig. [7^. When de- 
composing the AGN into host and nucleus, some part of the 
host's flux is systematically attributed to the central point source. 
However, this approach introduces fewer systematic effects than 
leaving the nucleus unaccounted for, as we will show below. The 
scatter in AC becomes larger toward lower H/N, independent of 
morphological disturbances. This explains why the C/A diagram 
of nucleus-subtracted AGN hosts is more noisy than that of qui- 
escent galaxies, and the effect is stronger for objects with lower 
H/N. 

The asymmetry index (Fig. fT2b ) is the only descriptor that 
shows a clear systematic effect: toward lower H/N, A s j m is in- 
creased, probably because of the artificial noise from the resid- 
uals after subtracting the nucleus. The scatter in AA (hence, 
the noise in the C/A diagram) also increases with decreasing 
H/N. This trend is much weaker for major mergers, however — 
mainly due to their high A m ; values — than for undisturbed 
galaxies. Since we are particularly interested in correctly identi- 
fying major mergers, this supports our results. 

The Gini coefficient (Fig. [TZb ) shows a trend similar to C, 
with slightly underestimated values of G at H/N > 0.5 and in- 
creasing noise toward lower H/N. The scatter is slightly smaller 
for major merger candidates. The M20 index (Fig.fT2h) is slightly 
overestimated for high H/N and underestimated for low H/N; 
the noise increases with decreasing host-to-nucleus contrast, in- 
dependent of the morphology. 

Fig. [O] demonstrates how the morphological descriptors are 
affected when an optical nucleus is not accounted for. Again, we 
computed the difference between the initial value (index ini) and 
that after the addition of a synthetic nucleus (index unsub) for 
C, A, G, and M20. Unsurprisingly, the concentration index (up- 
per left graph) becomes systematically higher due to the central 
point source, resulting in high negative values of AC U . This ef- 
fect becomes stronger toward weaker host-to-nucleus contrasts 
and affects major mergers stronger than undisturbed galaxies. 
The upper right plot shows the change in asymmetry index A. 
Particularly interesting here is the morphological dependence: 
while A on average becomes larger for undisturbed galaxies af- 



ter the addition of a nucleus (AA U < 0), the asymmetry of ma- 
jor mergers is slightly underestimated. Overestimated asymme- 
try indices probably are caused by nuclei that are not perfectly 
symmetrical but show a more complex shape of the PSF. The 
Gini coefficient shows the same trend as C, but without a notice- 
able dependence on morphology. M20 is systematically under- 
estimated (i.e. overestimated in absolute values), and this effect 
is slightly stronger for major mergers. The central point sources 
probably mask out bright structures that are off-center but still 
relatively close to the hosts' centers. 

To summarize, we find that an unsubtracted nucleus has a 
strong systematic effect on all four morphological parameters 
considered here. In contrast, subtracting the central point source 
systematically only affects the asymmetry index A. This implies 
that comparisons between quiescent galaxies and AGN with 
non-negligible nuclei will be biased if the impact of the nuclei is 
ignored. Similar results have been found e.g. by Sanchez et al. 
(2004) or Gabor et al. (2009). 

5.3. AGN host colors 

We constructed our quiescent galaxy sample in such a way that 
a fair comparison to the real AGN hosts is feasible: they are 
matched in their redshift and apparent magnitude distributions. 
A lower limit on the Sersic index was applied to reject bulgeless 
disks among the non-AGN. The initial selection criteria did not 
use any information on spectral type, except that galaxies clas- 
sified as AGN within the COMBO- 17 survey did not enter the 
comparison sample. 

For testing purposes, we used the observed color as an addi- 
tional constraint for the definition of the non-AGN sample. This 
analysis is restricted to the GEMS part of our sample, which of- 
fers v- and z-band HST/ACS imaging, while the STAGES survey 
only comprises v-band data from HST Note that the nucleus- 
subtracted HST images are our only option to derive host col- 
ors — the ground-based COMBO- 17 data have a totally insuffi- 
cient resolution for a host-nucleus decomposition. 

Our sample holds 12 real AGN hosts with available GEMS 
data. Their colors span the range 0.86 < (v — z) < 2.4. Matching 
the comparison sample to this color distribution reduces it to a 
total of 62 galaxies. 
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Fig. 12. Illustration of the impact of a subtracted nucleus on the quantitative morphologies of the simulated AGN galaxies in 
concentration index C (upper left), asymmetry A (upper right), Gini coefficient G (bottom left), and Moo (bottom right). As a 
function of H/N, all plots show the initial parameter value before (index ini) and after addition & subtraction of a synthetic nucleus 
(index sim), normalized to the initial value. See text for details. 



We repeated the KS tests using this new non-AGN sam- 
ple matched in redshifts, magnitudes, and colors. Comparing 
the AGN to major mergers, undisturbed, and unclear objects 
in C/A space, we now find P mm = 0.04, P unr ji s = 0.24 and 
^unci = 0.03, respectively. As in our initial analysis, undisturbed 
galaxies show the highest probability of stemming from the same 
parent distribution as the real AGN hosts. In G/M20 space, the 
respective probabilities are P mm = 0.29, P U ndis = 0.59 and 
^unci = 0.18. Again, undisturbed galaxies are most similar to the 
real AGN hosts. Compared to the initial results without a color 
selection, the difference in probability between undisturbed and 
major merger cases is reduced. 

We emphasize that, because of small number statistics, these 
results are not as robust as our initial analysis. The additional 
color criterion shrinks the comparison sample of major mergers 
from 15 objects to just six. In particular, four very blue major 
mergers with colors (v— z) < 0.9 are excluded from the quiescent 
sample. A visual inspection of these galaxies does not reveal ob- 
vious morphological differences to the rest of the major merger 
cases. But the very blue ones show relatively low absolute M20 
values and excluding them slightly affects the KS test results in 
G/M20 space (by increasing the value of P mm ). 

Nevertheless, the reanalysis with a color-constrained com- 
parison sample confirms our initial results. The fact that major 



mergers have much bluer average colors (((v - z)) = 0.95) than 
the real AGN hosts (((v — z)) = 1 -65) can be taken as some spec- 
trophotometric evidence that major mergers are scarce among 
the host galaxies. This agrees with previous findings by, e.g., 
Sanchez et al. (2004) or Schawinski et al. (2007): at least up to 
redshifts z ~ 1, AGN hosts show intermediate — but not very 
blue — colors. 



5.4. Discussion 

Our quantitative morphological analysis shows that the majority 
of 21 AGN in our sample reside in only mildly disturbed or even 
undisturbed galaxies. The same result is found in a visual clas- 
sification of the nucleus-subtracted images (see Table 2): Ten 
of them (48%) have an undisturbed morphology, 11 (52%) are 
unclear cases, but none has been classified as a major merger. 
For comparison, the quiescent sample comprises 15 (7%) major 
mergers, 95 (47%) undisturbed and 93 (46%) unclear objects. 

Hence, neither with a quantitative analysis nor by visual in- 
spection we find evidence that a significant fraction of the AGN 
hosts are undergoing major mergers. This is at variance with a 
naive expectation for a connection between major mergers and 
AGN. In principle, there are three possible explanations for this: 
1. Morphological descriptors like C/A or G/M^o might miss a 
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significant fraction of strong gravitational interactions on single- 
orbit ACS data of intermediate-redshift galaxies; 

2. there might be a delay between a major merger and the onset 
of the optical AGN phase: the merger signatures could fade and 
become undetectable (or very hard to detect) before the nuclear 
accretion starts, or 

3. violent gravitational interactions might not be the main trig- 
gering mechanism of nuclear activity in the redshift regime con- 
sidered here, which corresponds to look-back times between ~ 5 
and ~ 8 Gyr. 

Regarding the first of these possibilities, our analysis shows 
that major mergers should be detectable at a high success rate, 
even in the presence of residuals from a subtracted nucleus. This 
is what we conclude from the KS tests between undisturbed and 
strongly interacting galaxies (probability P = 8.6 x 1CT 4 using 
C/A or P = 1.1 X 1(T 3 using G/M 20 , cf. Fig.EJand TableD as 
well as from the discussion of the trends in Fig. [12] However, 
even major mergers are triggering the nuclear activity for a 
significant fraction of our sample, it is not necessarily the case 
that the morphological imprint of the merger is still detectable 
by the time the AGN phase has started. The start of the SMBH 
accretion and the time of the strongest morphological perturba- 
tions (hence the best "detectability" of the merger) might not 
occur simultaneously. 



This brings us to the second option listed above. 
Observational evidence for a delay between a merger and the 
onset of an AGN phase has been found e.g. by Schawinski et 
al. d2007b and Bennert et al. (120081 ) . Synthetic observations to 
estimate the time during which a major merger could be identi- 
fied have been carried out by Conselice (2006; using the CAS 
space) and Lotz et al. (2008; based on CAS and G/M%y). The 
involved time scales depend on very many parameters, such as 
the mass ratios, gas and dust content, or the geometry and sense 
of rotation of the merging galaxies. According to these authors, 
the phase of the strongest perturbation occurs at the earliest a 
few lOOMyr and at the latest > 1 Gyr after the first encounter. 
However, these analyses did not attempt to estimate the time 
when the nuclear gas accretion starts, nor did they include the 
impact of an AGN's optical nucleus on the image analysis. 

Using numerical simulations, Hopkins et al. (2008 ) have in- 
vestigated the various evolutionary stages of merger-induced nu- 
clear activity. The encounter and major merger phases are fol- 
lowed by the coalescence and start of the SMBH accretion. By 
the time the object would be observed as an AGN at rest-frame 
optical wavelengths, the tidal features might have faded. For dis- 
tant objects, the merger signatures could be only detectable with 
very deep observations, but not with single-orbit ACS data such 
as those from GEMS or STAGES. 
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Concerning the third hypothesis, it can be argued that at least 
at higher redshifts and luminosities, there is evidence both from 
observations and numerical simulations for a connection be- 
tween major mergers and AGN phases. Ultraluminous infrared 
galaxies (ULIRGs), e.g., which are mostly found in merging 
systems (e.g. Sanders & Mirabel 1996), have a much higher 
AGN fraction than the overall population of luminous galaxies 
(e.g. Canalizo & Stockton 2001). Simulations were used to show 
that major mergers can efficiently transport gas to the central 
parts of the merger remnant to feed the central SMBH (e.g. Di 
Matteo et al. 2005'). At later stages, the feedback by the AGN 
can terminate the star formation; this scenario might be a ne- 
cessity to produce the "red and dead" ellipticals observed in the 
present-day universe (e.g. Khalatyan et al. 2008). 

Studies using very large samples drawn from the Sloan 
Digital Sky Survey, on the other hand, did not find an increased 
fraction of interacting galaxies among AGN hosts at low z 
(e.g. Kauffmann et al. 120031 Li et al. 120081 ). It has to be noted 
that these local studies on average cover much lower AGN lu- 
minosities than intermediate-redshift samples like our own. A 
straightforward comparison is thus not feasible. The observa- 
tional results at differing redshifts are diverse but indicate that 
a possible connection between major mergers and nuclear ac- 
tivity seems to weaken toward lower redshifts and lower AGN 
luminosities. 

Our results are consistent with weak gravitational inter- 
actions (fly-bys, minor mergers) playing an important role in 
AGN triggering at z ~ 0.7. This would also agree with the in- 
creased galaxy density in the vicinity of quasars (e.g. Hennawi 
et al. 120061) and the enhanced AGN fraction in galaxy pairs 
at low (e.g. Ellison et al. 1201 II) and intermediate redshifts 
(e.g. Silverman et al. 1201 II ). Tidal interactions or minor merg- 
ers can induce the transport of gas with low angular momen- 
tum to the central parts of the host to (re-)start SMBH accretion 
(e.g. Hernquist & Mihos 1995); the host morphology would be 
much less affected than in a major merger event. Indeed we find 
that the AGN hosts are most similar to undisturbed galaxies, in 
the sense of highest KS test probability (P = 0.32 or P = 0.22, 
depending on which pair of descriptors is used). Our visually de- 
rived fraction of undisturbed AGN hosts is 48%, and it is likely 
that more objects with only weak tidal features are among the 
remaining 52% of unclear cases (Table[2]i. 

The lack of signatures for frequent strong interactions is 
consistent with the findings at similar redshifts by Grogin et 
al. (120531), P ierce et al. J2007b . Gabor et al. ( [20091 ), and Cisternas 
et al. d201 II ). One major difference between the first two of these 
and our analysis is that they did not account for the potential 
influence of an optical nucleus. The first three of the above sam- 
ples mainly comprise type-2 AGN since they are X-ray selected 
(Gabor et al. give a type-1 AGN fraction of ~ 10%). Only the 
data set by Cisternas et al. d201 II ) mainly consists of type-1 AGN 
(contributing 60 %); our own sample is made up entirely of type- 
1 AGN. Despite these differences in the type mix, our results 
and these previous studies agree qualitatively. This is to be ex- 
pected in the context of the classical unification scheme of AGN 
(e.g. Urry & Padovani l 19951 ), where the distinction between type 
1 and type 2 is assumed to be an observational effect and not an 
intrinsic property. 

To further interpret the observational findings, it would be 
highly desirable to create synthetic images based on numeri- 
cal simulations. With simulations that include galaxy interac- 
tions and mergers in various configurations and a prescription 
for AGN activity (as, e.g., Di Matteo et al. 120051 ). sets of syn- 
thetic observations would have to be generated to investigate, for 



a given data quality, the detectability of interaction signatures 
during the AGN phase. Such an approach is, however, beyond 
the scope of this paper. In the future, we will also aim at devel- 
oping a morphological descriptor scheme that is better capable 
of tracing weak interaction signatures than those we currently 
have at our disposal. 

6. Summary 

We carried out a quantitative comparison between the morpholo- 
gies of the host galaxies of AGN and quiescent galaxies at red- 
shifts 0.5 < z < 1.1. The imaging data were taken from the 
large HST/ACS mosaics of the GEMS and STAGES surveys. 
Our main aim was to test whether the AGN host galaxy popula- 
tion shows an increased fraction of morphological perturbations 
with respect to non-AGN. We created synthetic images to in- 
vestigate the impact of an optical nucleus on the morphological 
analysis of AGN host galaxies. To quantify galaxy morpholo- 
gies, we used the asymmetry index A, the concentration index 
C, the Gini coefficient G, and the M20 index. A sample of ~ 200 
synthetic AGN was matched to 21 real AGN in terms of redshift, 
host brightness and host-to-nucleus ratio to ensure a reliable 
comparison between active and quiescent galaxies. The optical 
nuclei of type-1 AGN strongly affect the quantitative morpho- 
logical properties of the underlying host galaxy. When these ef- 
fects are accounted for, active and inactive galaxies show similar 
distributions in C/A and G/M20 space. The morphologies of the 
AGN hosts are clearly distinct from galaxies undergoing strong 
gravitational interactions. This is at variance with a naive ex- 
pectation for a connection between major mergers and AGN+. 
However, the signatures of gravitational interactions might be 
too weak to be detected by the commonly used descriptors like 
C/A due to a delay between the strongest morphological pertur- 
bation and the onset of the AGN phase. 
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